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Several numerical examples are presented to show the predominant features of radio 
wave propagation over an inhomogeneous earth. These are based on the theories derived 
previously [Furutsu, 1957a, 1957b, 1959, 1963] in which the height and also the electrical 
properties of the earth's surface were assumed to change discontinuously several times along 
the wave path; thus, the terrain represented could include ridges, cliffs, bluffs at a coastline, 
etc. The theory is briefly reviewed, and numerical results are presented for the spherical 
earth approximation and are compared to those for the flat earth approximation. For a 
perfectly conducting flat earth, there are well-known formulas available in terms of the 
Fresnel integral, and the spherical earth results are compared to those obtained using 
these formulas to show the agreement at short distances. A few interesting phenomena are 
also illustrated, such as the obstacle gain due to a ridge on a lossy ground and the variation 
of field strength caused by a change of receiver (or transmitter) height when the wave is 
propagated over a mixed path. Finally, sets of graphs are also included to aid in evaluating 
the effects of a ridge or a cliff on a homogeneous earth ; they can be used when the propagation 
distances arc sufficiently large on each side of the ridge or cliff. 

1. Introduction 

The purpose of this paper is to determine the field strength when a radio wave is propa- 
gated over a terrain varying both in height and in electrical properties. 

The problems of mixed paths over a smooth earth have been investigated by many authors 
[Bremmer, 1954; Clemmow, 1950; 1953; Feinberg, 1944; 1946; 1959; Godzinski, 1958; Wail, 
1956; 1957; 1961], and one of the authors of this paper [Furutsu, 1955a; 1955b; 1955c]. Some 
of their results were given in terms of a convolution integral which is usually evaluated by a 
numerical method. Mixed paths mean that the wave path is over ground consisting of several 
sections of different electrical properties, but in these papers, no change of height of terrain 
had been taken into account. 

However, it had been the main purpose of some papers by Furutsu [1957a; 1957b; 1959; 
1963] to determine the effects on the groundwave propagation over an inhomogeneous earth 1 
caused by significant changes in the height of terrain, such as ridges, cliffs, etc., and also how 
these effects could be treated. The model of the terrain, the height of a spherical earth surface 
and its electrical properties were assumed to change discontinuously along the wave path 
several times. This is illustrated in figure l. 2 The formula of field strength for this model was 
obtained in the form of a multiple residue series. This formula reduces to that for mixed paths 
in the special case of a smooth earth and further to the ordinary Van der Pol-Bremmer series 
in the case of a homogeneous earth. 

The basic terrain used in this paper is illustrated in figure 2. This figure shows the terrain 
consisting of two sections of different earth radii, a 2 and a 4 , and different propagation constants, 



1 Sec, however, papers by HufTord [1952] and also by Wait [1963]. 

2 In this figure, a„is the radius of the nth section, r n is its length, i\ and x n +i indicate the transmitter and receiver locations. 
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Figure 1. The form of the terrain used in the theorx 
of inhomogeneous earth. 




Figure 2. The form of the terrain consisting of two 
sections for (3). 



k 2 and & 4 , respectively, with a ridge of radial distance a 3 (measured from the earth's center to 
the top of the ridge) at the boundary of the two sections. The transmitter and the receiver 
are on each side of the ridge at the points Xi and x 5 , whose radial distances are z x and z 5 , respec- 
tively. Their distances from the ridge are r 2 and r 4 , being measured along a mean earth surface 
of radius a. 

The attenuation coefficient A is denned by (A.2.1), 3 i.e., if E is the field strength to be 
obtained, 



E=2AEn, 



E - 



47r(/' 2 +r 4 ) 



e-ttiOvW, ki=27r/\ : 



where X is the wavelength in free space and E may be regarded as the field strength in free 
space. The attenuation coefficient A thus defined is conveniently given in terms of the numeri- 
cal distances c 2 and c 4 defined by 



c 2 = (r 2 /a) (ha/2) 1 *, c,= (rja) (ha/2)^, 
and also the numerical heights y u y 2 , y 3 , y 4 , and y 5 defined by 

y 1 =h(z 1 -a){2/k 1 a) 1 ^, y 2 =h(a 2 -a)(2/k 1 a) 1 ^ ) y^=k 1 (a 3 -a)(2/k 1 a) 1 ^ J 

y,=k 1 (aA-a)(2/k 1 a) 1 ^\ y 5 =k 1 (z 5 -a)(2/k 1 2) 1 ^. 
Thus, using (A.3.1), A can be given in the form 

A=^2 { (c 4 +c 2 )/c 4 } ll2 A{y^ c 4 ) h T™ fe) %4f t J t2 {y l2 ) , 
with the conditions 

Oi— a 2 )/r 2 «l, (a 3 — a 2 )/r 2 «l, (s 5 — a 4 )/r 4 «l, (cx 3 — a 4 )/r 4 «l, 

Vij=yi-yj=-yji, 

A(y^c A ) t ^{TTc A ) ll \t,-(il)- l j t ^y bi )e^ [— i{c 4 (»4+0+'/4}], 



where 



a) 



(2) 



(3a) 



(3b) 

(4) 
(5) 



T (3) (c 2 ), 4 . (2 = {qJ' u (y Si )fi 2 (y32)-zJ ti (yu)f, 2 (yz2)} {2/4— 2/2+^4—^2} _1 (*»— jST'expl— ifc(iM-*i)]- 



Here, the set of values t m (m=2, 4) stands for the roots of the equation 

W / '(U-2m^(f m ) = 0; 
W'(i) is the first derivative of the function W(t) defined by 

W(-t) = (irt/3) 1/2 e- i2 *< 3 m%(it 3 < 2 ), 



(6) 
(7) 
(8) 



3 Equation numbers preceded by an A refer to equations in Part I of Furutsu [1963]. 
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and 

/7 , N1/ , J h-y/k 2 rn -k 2 Jk 2 m , Vertical Polarization 

^=(&ia/2) 1/3 X-< > > W 

I -yjk 2 m —k 2 jk h Horizontal Polarization. 

The function f tm (y) is the ordinary height-gain function defined by 

J tm (y) = W(t m -y)/W(t m ), (10) 

and the function j\ n {y) is defined by 

iJ' , m (y) = - W*yy, m (y) = W'(t m - y )/W(t m ). (ii) 

Thus, it can be seen from (7) that 

/'. m (0)=/U0) = l. (12) 

In the special case where az=a 2 {a 2 >a A ) or a 3 =a i (a 2 <a i ) , the terrain would represent 
a cliff, and then 

f 24/! 4 (2/24) —^4(2/24), 2/2>2/4 

Xexp[-ic 2 (2/ 2 +^ 2 )]X^ ... U3) 

(^ 247*2(2/42) —g 2 7* 2 (2/42) , ?/ 2 < 2/4. 

In the case of a 2 =<x 4 and q 2 =q^ the terrain would represent a ridge on a homogeneous 
ground. In this case, both the numerator and the denominator of 

vanish for £ 4 =/ 2 ; then 

T< 3 > (o 2 ) «,, „= [ { 1 - y 32 fe-al)" 1 }/f 2 (2/ 32 ) 

(14) 

+^fe-(/?)- 1 {y^(l/3 2 ) -/«"(?/:«)}] exp [-tcfo+tf]. 

Further, in the case of smooth earth where 2/2=2/3 = 2/4 = 0, (13) gives, on using (12) and (14) : 

T^(c 2 ) ti! , 2 =(q i -q 2 )(t i -t 2 )-\t 2 -ql)- i e- ic 2'2 (q i9 *q 2 ) 

' e -ic 2 H ti=t 2 

(<Z4=2 2 ) (15) 

^0, t 4 3**„ 

and thus, formula (3) becomes: 

A=^{*(c 2 +c i )} 1 Hq i -q 2 )(U-t 2 )- 1 (U-qt)- l (h-q 2 i )- 1 

X/ !2 (2/i)A(2/5) exp [— iicitz+CiU+w/i}], fe^*) (16) 
=S {xfe+c*)} 1 ' 2 ^-^) -«y, 2 (2/0^(2/5) exp[-t{(«,+C4)«.+T/4}], (24=22). (17) 

'2 
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Equation (16) is a special case of the established formula for mixed paths over a smooth 
spherical earth consisting of several different sections [Furutsu, 1955c], and (17) is the ordinary 
formula for a homogeneous earth. 

In the Van der Pol-Bremmer series for a homogeneous, smooth earth, the convergence 
becomes poor when the propagation distance is sufficiently short, and likewise, the convergence 
of the series (3) becomes very poor when the propagation distance on one or both sides of the 
ridge is very short. However, in the latter case, the effect of the earth's curvature is so small 
that the flat earth approximation can be used; the result in the flat earth approximation then 
provides the asymptotic form of the series (3) as a— >oo. 

In order to present the flat earth formula corresponding to (3), it is convenient to introduce 
the parameters defined by: 

d 2 =iqic 2 = — i(kilk£ 2 kir 2 l2, 

d A =iqlc i = — i(ki/k , A ) 2 k i rJ2, 

/i = K 1/2 2/i2^' W4 =(^i-« 2 )(/:i/2r 2 ) 1 / 2 6- /4 , 

/5=K 1/2 2/5 4 6- /4 =(^ 5 -a 4 )(i 1 /2r 4 ) 1/2 6^ /4 , 

f2=^ l/2 y 32 e^ / "=(a 3 -a 2 )(k l /2r 2 y /2 e i - / \ 

/4=K" 1/ V84e* W4 =(a»-a4)(*i/2r,) 1/2 e*' /4 , 

n 2 =r 2 /(r 2 +r 4 ), 

n 4 =r 4 /(r 2 +r 4 ), (18) 

with 



7 , J kll^kl-kl Ver. Pol. / 

k' m =l (m=2,4) (19) 

{kll^kl-kl Horiz. Pol. 

all of which are independent of the earth radius a. It may also be noted that d 2 and d 4 are 
the Sommerfeld numerical distances on the respective sides of the ridge. 

In the special case where both the transmitter and the receiver are on the ground, i.e., 
when Zi = a 2 , Zs=a A orji=f 5 =0, the attenuation coefficient A is given, with the aid of (A. 4. 2) 

to (A.4.4), by 

A^F{d 4 ,j 4 \d 2 J 2 )^e~^ + ^[g{^ 

X \^d 4 d 2 \n 4 n& (ft+i^di) g (/ 2 + W5) + (d 4 /n 4 ){8 (<r 4 , Un 2 j\+-yjn 4 f 2 )/a,) 
-<S(p 4 , (f*+ijdd/pd}+(d2/n 2 ){g(cr 2 , (^4+VM/<r2)-£(/>2, U+WS)//*)}]]. (20) 

Here the vertical line in the arguments of the function F defined above is used to distinguish 
between the variables on each side of the boundary, and 



(TA=^n A J 4 r-^n 2 J2+i^d A ln A , P4=i\d 4 n 2 /n 4 —f 2} 

<T 2 =^n 2 f 2 —-yJn 4 f A +i-\/d 2 /n 2 , p 2 =i-\Jd 2 n 4 /n 2 —f i , (21) 
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and 

g( 2 )=JLe* 2 J" e-* 2 dx, 

S(z,n/z)=(^)\" 2+ " 2 f dxe-' 2 C e~" 2 dy. (22) 

The analytic expansions for S (z) and £ (z, n /z)&re treated in the appendix of Furutsu [1963]. 
In the general case of Zi7 £ a 2 and z^a M 

A=F(d M MMd 2 , /2+/1) + %{F{d i} /4+y s |0, f a -fd -F(d t , / 4 +/ 6 | 0, / 2 +/i) 

+ ^(0, /*-/ 5 |4 J2+j\)-F(0, f 4 +Md2,f2+J l )}+%{F(0J i +j 5 \0, f i +fi+F(0,f t -j i \0,f a -fd 

-f(q, / 4 +/ 6 |o, / 2 -/i)-F(o, y 4 -/ 6 |o, /,+/i) } • (23) 

When the earth's surface is perfectly conducting and the wave is vertically polarized, (18) 
and (19) yield c/ 2 =rf 4 = 0. Thus, (20) becomes simply 

A = e ~(fl+fi)g UnJi+^nJ 2 ), 

v^+v^/ 2 ={^ 4 +«} y/^+k) eMi > (Zi=a » z *= a ^ (24) 

As is seen from the definition of the function & (z) in (22), (24) is given by the Fresnel integral 
and, as expected, is exactly the same as the result for ordinary knife-edge diffraction. 

On the other hand, when the earth's surface is smooth and both the transmitter and 
receiver are on the ground so that 

/l=/ 2 =/ 3 =/4=/5=0, 

formula (20) becomes 

-4=1— i V 7T ( -y/(U/n 4 + -vW^) ~ [-y/diddn^S (i-\/d 4 ) 8{i^d 2 ) 

+ (dt/n A ) { <§ (i^djrii) — 8 (i-yJd A n 2 /n 4) -y/njn 2 ) } 

+ G/2M2) {£ (i^da/fh) ~ 8 (i^d 2 njn 2 , ^n 2 /n A ) }]. (25) 

This equation is essentially equivalent to the formula first obtained by Clemmow [1953] and 
also the same as by Furutsu [1955a]. Further, if the earth is homogeneous so that (/ 2 /fh = d,Jn 4 , 
(25) becomes 

A=l-i(wd 2 /7i 2 y /2 6 (i^d^ 2 ) = l-i2(d 2 ln 2 ) 1/2 e-^n 2 f _ e — dx, (26) 

which is the formula commonly used for a homogeneous flat earth. In deriving this formula, 
the following relation is used: 

<S (a, b/a) + 8 (b, alb) = 8 (a) 8 (b). (27) 

In formula (3), it frequently happens that the propagation distance on one side of the 
ridge, say r 2 , is small enough so that c 2 <C < 0> while the propagation distance on the other side 
is large enough so that c 4 >\. Under these conditions the flat earth approximation may be 
used on the c 2 side, while keeping all spherical properties as they are on the c 4 side; i.e., when 

c 2 «l, |c 2 (y 4 2+*4)|«l, <Wia«l, c 2 t/ 32 «1, (28) 
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it has been shown [Furutsu, 1957b] that 

S^ 131 (<*).„ ,J t2 (yi2)^B ti (y 3i , d 2 , f 2 , /0 s[ { (iJqAf tl <3ti-f, t (3hd}Jdj(-&, /1+/1) 

+ K(Lj<L 2 )f u (yz i )Jd 2 {J(0, / 2 -/i)-J(0,/ 2 +/i) } +l/, 4 (2/34) { J(j2-ji)+J(f 3 +tt } ]«-'««■*»-»*.'., 

(29) 

where /i,/ 2 , and c? 2 are the same as defined in (18) and 

J(xJ)=e-* 2 {g(j)-g(f+ix)}/x, 

JU)=e- f2 gU)> (30) 

Furthermore, when 2/ 3 4<<l, 

/« 4 (2/34)^1— 2/3424=l+i(*l/*i)^l (fla — ^), ?/34<<l 

and 

(^/^^^fe/^)-^/^^ 2/34«i. (31) 

Also, when y 34 =0, the right-hand side of (29), B, (0, d 2 ,f 2 ,j'i), turns out to be independent of 
£ 4 and thus applying formula (29) to (3) on the condition (28), we get 

A~S^(lfc4, C4), 4 Bi 4 (0, *,/2,/i)=A| Ca . Bi 4 (0, dtjiji), ^|, 2 =o=2>%64, c 4 ), 4 . (32) 

*4 *4 

The latter series is the ordinary Van der Pol-Bremmer residue series for a homogeneous earth. 
In this paper, the above formulas are used to find the attenuation for several typical cases of 
inhomogeneous earth. Emphasis is first placed on the comparison of the numerical values in 
the spherical earth approximation with those in the flat earth approximation, since they are 
expected to asymptotically approach each other at short distances. Although there is almost no other 
literature available for wave propagation over a terrain of finite conductivity, there is the well- 
known formula given by the Fresnel integral for a ridge on a perfectly conducting plane, which 
therefore gives an interesting means of checking the corresponding spherical earth values. 
Secondly, a few interesting physical phenomena (such as the "obstacle gain" and a complicated 
variation of the field strength versus the height over mixed paths, etc.) are discussed with nu- 
merical illustrations. Finally, a set of graphs is presented for use in calculating the effects of a 
ridge and a cliff on the groundwave propagation; this is an extension of the work shown in 
figures 8 and 9 in [Furutsu, 1963]. They can be used in the same way as the ordinary height 
gain function for the transmitter and/or the receiver, when the propagation distances on both 
sides of the ridge are large enough so that 

c 2 -yU 2 -yii 2 »i, e 4 -2/ 3 f-^ 2 »l. (33) 

However, the above condition may be too strict for practical use, as will be seen by a few 
examples shown in the appendix. 

In the last section, the more general case is treated in which the inhomogeneous earth 
consists of several different sections with a ridge at each section boundary. This is illustrated 
in figure 12 and the formula for the attenuation coefficient is given, which was derived in earlier 
work [Furutsu, 1957b]; this formula corresponds to (3) for two sections with a ridge between 
them. 

2. A Ridge on a Homogeneous Earth 

The basic formula for the attenuation coefficient is given by (3) in the case of a spherical 
earth and by (20) or (23) in the case of a flat earth; for a ridge on a homogeneous earth, the 
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Figure 3. A ridge on a homogeneous earth {flat or 
spherical) for figures 4 o,nd 5. 
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0.02 0.05 0.1 0.2 0.5 1.0 

(r/a)(k,a/2)i 

l'"i<;ri;K 1. Comparison of the spherical and the flat 
earth values of attenuation due to diffraction by a 
ridge on a perfectly conducting earth (vertical 
polarization) . 

Attenuation versus the numerical distance (r/a)(/cia/2)i/3 when 
r 2 =r4=r, for several values of the numerieal ridge height y = k\h 
(2/Aio)V8. 



attenuation is found when yz2=Vu : =y = k\h{2lk l a) l,z and q 2 = q4 or k 2 = k±. For the numerical 
illustration the transmitter and receiver were placed on the ground at equal distances from the 
ridges, as shown in figure 3, and thus r 2 =r±=r or c 2 = c 4 =c=(r/a) (kia/2) 1/Z and Zi=-a 2 =a±=z b or 
2/12=2/54=0. 

2.1. Case of a Perfectly Conducting Earth 

When the earth is perfectly conducting and the wave is vertically polarized, q 2 = q 4 =0 
according to (9). Furthermore, if the earth is assumed to be flat, the attenuation coefficient A 
is given by the Fresnel integral of (24). In figure 4, the absolute value of A is plotted versus 

c=(r/a)(k 1 al2) lf * 

for several values of y and is displayed by the broken lines. 4 The corresponding values in the 
case of a spherical earth, shown in the same figure by solid lines, were obtained by evaluating 
formula (3) for g 2 =^4= — ^10~ 6 (which is effectively zero). As expected, the figure clearly 
shows that the spherical earth values asymptotically approach the corresponding flat earth 
values when c<<0- It may also be noted that the attenuation increases with the height of 
the ridge in the range of variables presented. 

4 Although both the variables y and c are dependent on the earth's radius a, the attenuation coefficient A is independent of a, as is seen from 
relation (18). 
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FLAT EARTH VALUES 

—o— SPHERICAL EARTH VALUES 



Figure 5. Comparison of the spherical and the flat 
earth values of attenuation due to diffraction by a 
' on an earth with a finite propagation constant. 
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Figure 6. Mountain gain and diffraction loss by a 
ridge on a lossy surface. 



2.2. Case of a Finite Conducting Earth 

The result of the comparison using the same data except for g_2=g_4= —i50 is presented in 
figure 5. As in figure 4, the spherical earth values approach the flat earth values when c<C<Cl. 
However, there is a remarkable difference between the characteristics of the two figures; i.e., 
in figure 5, the higher the ridge, the smaller the attenuation. This implies that the propagation 
loss along the lossy smooth surface is much greater than the diffraction loss by a ridge on the 
surface. Thus, we have the so-called obstacle gain. However, since the diffraction loss in- 
creases with a decrease of distance from the ridge, the obstacle gain is apparent only in a 
limited range. This situation is illustrated in figure 6, which is a continuation of the flat earth 
values in figure 5. 

On the other hand, when the propagation distances on both sides of the ridge are large 
enough to satisfy the conditions given in (33), the convergence of the series in formula (3) is 
fast enough so that only the first term of the series is necessary to give a sufficiently accurate 
value. This is reflected in figure 5 by the fact that the three solid lines for the spherical earth 
values are almost parallel to each other for the large distances occurring when c 2 =c 4 ^l. 
Therefore, in this case, the eifect of the ridge (of height h), located on a homogeneous spherical 
earth of earth constant g 2 , can be obtained by multiplying the attenuation coefficient for the 
homogeneous earth by the ridge gain factor T R (p) given by: 



T B {p)={l-y(t 2 -q%- 1 }My)+ql(tl-®- 1 {fAy)-f l t 



y=k 1 h(2/k 1 a) 1/ 



(34)' 



5 This equation is the same as (A.3.3) except that the second term is omitted in the latter. However, all the values in fig. 5 [Furutsu, 1959] 
and fig. 8 [Furutsu, 1963], were computed using the correct equation (compare with fig. 13 in this paper). 
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r 2 -i q 2 



Figure 7. A cliff or a bluff on an earth {fiat or 
spherical) for figures 8, 9, and 10. 



This is the term in square brackets on the right side of (14) where f 2 stands for the first value of 
t 2 for the homogeneous earth specified by q 2 \ for later convenience, p is defined as 

p=2- 1/s y=k 1 h(ha)- 1/ \ A=height of ridge, (35) 

and q 2 can be expressed as 

q 2 =2- m K- 1 e- i{b/2 ^ /4 \ (36) 

where if and b are the parameters introduced by Norton [1941]. 

A set of graphs is presented in figure 13 for \T R (p)\ versus p for various values of K and b. 

3. A Cliff or a Bluff at a Coastline 

When the obstacle is a cliff or a bluff (of height h), the basic formula in this instance is 
again given by (3) for the spherical earth and by (20) or (2:]) for the flat earth with 7/ 34 = and 
y 32 = kih(2/k i ay / '\ Further, if c 4 is chosen sufficiently large to satisfy the condition (28) for 
t A , formula (29) is available for the range of c 2 <C < 0« The latter is the case of figures 8b, 9, 
and 10 to be explained in the following, where the absolute values of the ratio A/(A\ C2={) ) are 
plotted versus c 2 for several values of the parameters. 

For the numerical illustration, the antennas were again placed on the ground (zi=a 2) 
z 5 =a A =a s , or y 12 =y 5i =0) , as shown in figure 7. 

3.1. A Cliff 

In this case, k 2 =k 4 or q 2 =q 4 and, when the earth's surface is perfectly conducting and the 
wave is vertically polarized, q 2 = q 4 = 0. Then the attenuation coefficient is again given by the 
Fresnel integral of (24) with/ 4 = for the flat earth approximation. These values are shown 
in figure 8a and 8b by the broken lines, and the corresponding spherical earth values are shown 
by the solid lines (as in section 2, q 2 =q± was again chosen to have a very small value — i!0~ 6 , 
which is effectively zero). Figure 8a shows the absolute values of A versus (r/a) (&ia/2) 1/3 , when 
r 2 =r i =r, while figure 8b shows the relative values of \A\ versus the numerical distance 
c 2 =(r 2 /a)(k 1 a/2) 1/3 , when c 4 = l0; the latter values of |^4| were computed using formula (32). 

A similar comparison is presented in figure 9 for the same data except for q 2 = q 4 = —i5Q, 
where the flat earth values were computed using the function B t4 (y u , d 2 , f 2 , fi) defined in (29) 
with 2/34 =/i=0 and q 2 =q 4 . It may be noticed here that, when z x = a 2 orf 1 = 0, B H is independent 
of the value of q 2 =q 4 on account of the identity (12). Thus, the flat earth values are the 
same as in figure 8b. 

As in the case of a ridge on a homogeneous earth, when the propagation distances on both 
sides of the cliff are large enough to satisfy the conditions in (33) for y u =0, the effect of a cliff 
(of height h) can be obtained by multiplying the attenuation coefficient by the cliff gain factor 
T c (p). On referring to (13) and figure 7, the factor T c (p) is found to be 

Tc(p)=i2(t° 2 -qt)- 1 y- 1 {f, l (y)-P 2 (y)}, (37) 

where p=2- 1/3 y=k 1 h(k l a)- 1/s . 
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Figure 8. Comparison of the spherical and the flat 
earth values of attenuation due to diffraction by a 
cliff on a perfectly conducting earth (vertical polari- 
zation) . 

(a) Attenuation coefficient versus the numerical distance (r/a) 
(/cia/2)i/3 when r2=ri=r. 

(b) Relative values of the attenuation coefficient versus C2= 
(r 2 /a)a-ia/2)i/3 when t 4 =10. 
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Figure 9. Comparison of the spherical and the flat 
earth values of attenuation due to diffraction by a 
cliff on an earth with a finite propagation constant. 

Relative values of attenuation versus the numerical distance 
(n!a)(hal2)U* when c 4 = 10, for several values of the numerical cliff 
height y=*kih(2/kia)W. 
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Figure 10. Comparison of the spherical and the flat 
earth values of attenuation due to diffraction by a 
bluff on a coastline. 



Upon reference to (31), we see that T c (p) tends to 1 as p— >0. 

A set of graphs (fig. 14) is included for \T c (p) | versus p for various values of K and b as 
defined by (36). These figures together with figure 13 for \T R (p)\ may be used to determine 
the gain (or loss) by a cliff or a ridge on a homogeneous earth; they are used in the same way 
as the ordinary height gain function to evaluate the gain caused by the variation in height of 
the transmitter and/or the receiver. 



3.2. A Bluff at Coastline 

This is the case where I^K^I&I ^ or vertical polarization or |</2|>>|<Z4| for horizontal 
polarization (see fig. 7). For the numerical illustration, q 2 and # 4 were chosen to be q 4 = — i50 
and q 2 =— H0~ 6 and the values obtained are shown in figure 10. The flat earth values were 
computed using formula (32) with 



y*2=y=kM2/k 1 a) 



1/3 



h being the height of the bluff. Figure 10 again shows that the flat earth values give the 
asymptotic values for c 2 <C < Cl. Also shown here is the interesting phenomenon of increasing 
field strength with the increase in distances from the bluff on the q 2 side (sea). 
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Figure 11a. The form of terrain used in conjunction 
with figure lib. 
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Figure lib. Field distribution over mixed paths 
consisting of two different sections. 

Attenuation versus the numerieal height kih{2lha) l ^ol the receiver 
for several values of the numerical distance C2. 



4. Field Distribution as a Function of Height Over Mixed Paths 

The field distribution as a function of the height over a homogeneous earth has been 
numerically treated by many people. However, there have been h*w attempts made to find 
the field in the ease of mixed paths. A sketch of mixed paths is shown in figure 1 la. In this 
sketch, the transmitter and the receiver are located over sections of different electrical properties 
and one of them, say the receiver, is allowed to vary in height. 

In figure lib, the field strength is plotted versus 6 



for c 4 = 10 and some values of 



k 1 h(2/k 1 a) i 

c 2 =(r 2 /a)(k 1 a/2) }i 



both for this case (solid lines) and for a homogeneous earth (dotted lines) specified by </.\. When 
the receiver is above the optical boundary point (noted in both figs. 1 la and lib by a cross), 
the field strength becomes more and more dependent on the earth's electrical properties on the 
transmitter (r/ 4 ) side; this is indicated by the asymptotic behavior of the curves in figure lib. 
It may be remarked that, although the field strength approaches the values for homogeneous 
earth at large heights, the points of approach are much higher than the optical boundary points; 
furthermore, the field strength oscillates with the height at larger heights. This oscillation 
may be reasonably interpreted as the interference between the principal wave and the wave 
induced at the boundary of discontinuity between the two sections of the inhomogeneous earth. 

5. Formula for the More General Case 

Formula (3) gives the attenuation when the wave is propagated over an inhomogeneous 
earth consisting of two different sections. This formula can be extended for the more general 
case, illustrated in figure 12. The terrain in this figure consists of several sections (each of 
which is a homogeneous surface having different heights and electrical properties), with a ridge 
at each section boundary [Furutsu, 1957b]. The transmitter and the receiver are located at 
the points X\ and x n+ i, and the attenuation coefficient A is then given by: 

A = ^t n t n _ 2 t t if t 2 (c/c n ) /2 A(y n+ltn , c^) tn T {n - l) {c n ^ tnt tn _ 2 

XT^\c n ^) tn _ 2ttn _ A . . . T^(c,)^ h T^(c 2 ) tit J t2 (y 12 ) (38) 



3 h denotes the height of the receiver. 
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Q n-I X n+ , 




P'igure 12. The general form of the terrain as used 
for (38). 

with the conditions similar to (3b). Here, the functions 

Myn+i,n,c n ) tn and T {m - l) {c m - 2 )t m ,t m ._ 2 

(m=n, n—1, . . ., 4) and ft 2 (yi2) are the same as those denned in (5), (6), and (10), 
respectively, and c is the total numerical distance. 

The above formula, (38), agrees with (3) in the special case of n = 4. Also, when a n = 
a n - l = a n -2= . . . =a 3 =a 2 = a, formula (38) reduces to the established formula for mixed paths 
on a smooth earth [Furutsu, 1955c]. Further, the following lemmas hold: 



lim S^W^/^fe 



a)t n ,t q — *S 



T (m \c g ) t 



a m ^a p 



lim J^ r>\c q ) tn , t J t (y sq ) = f tn (y sn ), 2 s ^a p . 



a m ^a p , 



C fl -» + tq 



(39) 



(40) 



These are required from the self -consistency of the formula when the length of any section 
involved in the inhomogeneous earth approaches zero, or when the receiver (or the transmitter) 
approaches the boundary of the section. 



6. Summary and Discussion 

Some numerical results, based on theory developed in previous paper [Furutsu, 1957a, 
1957b, 1959, 1963], are presented for several typical cases of inhomogeneous earth. The 
results are found for both the spherical earth approximation and the flat earth approximation, 
and, as expected, the spherical earth values approach the flat earth values at short distances. 
In the special case of a ridge on a perfectly conducting flat earth, the numerical results are also 
favorably compared with those obtained using the Fresnel integral. 

Many problems which seem to be independent (e.g., mixed paths on a smooth earth, 
terrains including ridges, cliffs, and bluffs at a coastline, etc.), can be treated by special 
applications of the same formula. For the first illustration, the effect of a ridge on a lossy 
homogeneous earth is presented, which shows the obstacle gain at large distances, and the 
diffraction loss at small distances. Using the same formula in the second illustration, the 
field strength is found for a varying receiver (or transmitter) height when the wave is prop- 
agated over a mixed path. These values are presented along with the values for a homogeneous 
earth, and the comparison is quite interesting, especially from the optical viewpoint; they 
approach the optical values at much higher points than expected and also show an interference 
with another wave which is supposedly induced at the boundary of discontinuity of the mixed 
path. 

The effect of locating a ridge or a cliff on a homogeneous earth is also found by the same 
basic formula, and sets of graphs (figs. 13 and 14) are presented showing this effect when the 
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Figures 13. Ridge gain function |T R (p) | as a func- 
tion of the numerical height p. 
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Figure 13. — Continued. Ridge gain function |T R (p) | 
as a function of the numerical height p. 
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Figure 13. — Continued. Ridge gain function |Tr(p)| 
as a function of the numerical height p. 
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Figure 14. Cliff gain function |T c (p)| as a junction 
of the numerical height p. 
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Figure 14. — Continued. Cliff gain function |T c (p)| 
as a function of the numerical height p. 
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Figure 14. — Continued. Cliff gain function |Tc(p)| 
as a function of the numerical height p. 



distance on both sides of the obstacle is large. As the graphs indicate, the field strength 
increases continuously with height (denoted by p) when the earth is a poor conductor (K«l), 
but decreases for small ridge (or cliff) heights as the earth becomes a better conductor. It may 
also be noted that for very large heights the gain caused by a ridge is much larger than that 
for a cliff of the same height. These sets of graphs may be used in the same way as the 
ordinary height gain function has been used for evaluating the effect of antenna height of the 
receiver and/or the transmitter. 

A subsequent paper will be published containing sets of graphs to aid in finding the field 
strength for various values of earth constants, ridge or cliff heights, propagation distances, etc. 



The authors express their thanks to L. A. Berry for the use of a few of his subroutines for a 
high-speed computer, and also to Mrs. J. E. Herman and W. S. Stevenson for their valuable 
assistance in obtaining the numerical results in this paper. 



7. Appendix 

In this appendix a few examples are presented showing the use of the set of graphs (fig. 
13) for the ridge gain function. 

The basic variables are the numerical distance c, the numerical height y=2 1/3 p, and the 
earth constants K and b defined in (1), (2), (35), and (36), respectively. In terms of more 
practical variables, they are given as follows: 

c=2.188X/ 1/3 a- 2/ V, 

2/=9.579X10- 3 X/ 2/3 a- 1/3 />, 

p=2- 1/3 2/=7.603X10- 3 X/ 2/3 a- 1/3 A, 
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e-1) 2 +s 2 }" 1/4 Ver. Pol. 
K=0.3627X(/a 

- 1/4 Hor. Pol. 



b = 



f(e 2 + . 2 ) 1/2 {(e-i: 

{{{e-iy+s 2 }-^ 

{2 tan" 1 («/*)— tan" 1 ! (e— 1)/«} Ver. Pol. 
(A.l) 
L80°- tan" 1 {(«—!)/*} Hor. Pol., 



with 



«=18<r//. 



Here / is the frequency in megahertz, a is the effective radius of the earth in kilometers, r is 
the propagation distance in kilometers, h is the height in meters, e is the relative dielectric 
constant of the earth (referred to a vacuum), and a is the conductivity of the earth in mil- 
limhos/meter. 

In the following, referring to figure 2, a ridge (of height 100, 200, or 300 m) is located on 
a homogeneous earth of effective radius 8500 km having the constants e=10 and <r=0.1 (10~ 15 
emu) and a vertically polarized radio wave of frequency 300 MHz is propagated from the 
transmitter to the receiver, both having the height 10 m, over the distances r 2 =r 4 =100 km. 
Then, by using (A.l), the following values are obtained for the necessary parameters: 

i£=0.00885, 

6 = 89.99°, 
2/12=2/54=0.210, 
2/32=2/34=2.10 (100), 4.21 (200), 6.31 (300), 

p=2- 1/3 2/ 32 =1.67 (100), 3.34 (200), 5.00 (300). (A.2) 

Here the values in parentheses are the corresponding ridge heights in meters. 

The exact values of attenuation for the data in (A.2) are found, by using formula (3), to 
be 

|^4|=6.20X10- 7 (100), 4.27X10" 6 (200), 1.91XHT 5 (300), (A.3) 

where the values in parentheses are again the corresponding ridge heights in meters. 

On the other hand, when there is no ridge, the corresponding attenuations can be obtained 
by using the Bremmer series for a homogeneous earth, and the method for obtaining these 
homogeneous earth values lias been given [Norton, 1941]. The homogeneous earth value, 
\A \, is then found to be 

|A | = 1.49X10- 7 . (A.4) 

Thus, the ridge gains |-4/^4 | are 

\A/A \=4:.l7 (100), 28.7 (200), 128. (300). (A.5) 

Here, if B is defined by 

B^c-fli-flh (A.6) 

giving the left side of condition (33), it has the following values for the data in (A.2) : 

5=1.61 (100), 1.01 (200), 0.55 (300). (A.7) 

Since &==90° according to (A.2), the ridge gains may be obtained from figure 13.6, and 
the values in (A.5) are found to be the same as those from the graph within the range of graphical 
error. Indeed, the exact values are 

4.22 (100), 30.2 (200), 147 (300). (A.8) 
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Here it is noticed that, as is seen from the values of B in (A. 7), the condition given by (33) is 
not satisfied for the above data used. However the values in (A. 8) may be sufficiently accurate 
for practical purposes. 

It will be worth while to check a case where B is negative; e.g., when r 2 =r A =50 km or 
c 2 =c 4 = 1.76 and the other parameters are the same as in (A. 2), the exact value of the ridge 
gain for h=100 m is found to be 

|A/A | = 3.52 (100) (A.9) 

with 

£=-0.5. 

Here, the corresponding ridge gain from figure 13.6 is about 4.4, and therefore the error caused 
by the use of the graph turns out to be about 25 percent or 2 dB. 
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